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A generalized Dirac equation is derived in order to describe charge carriers moving in corrugated 
graphene, which is the case for temperatures above 10° K due to the presence of flexural phonons. 
Such interaction is taken into account by considering an induced metric, in the same spirit as the 
general relativity approach for the description of fermionic particle moving in a curved space-time. 
The resulting equation allows to include in a natural way the presence of other phonon branches as 
well as an external electromagnetic field. It also predicts non-linear effects which are not present 
in the usual vector potential approximation used in most of publications on the subject, as well as 
the possibility of controlling electronic conductivity using pure sinusoidal strain fields. The non- 
linear terms are important at high temperatures, and can also lead to interesting effects, like e.g. 
resonances between flexural phonons and external electromagnetic fields. 

PACS numbers: 



I. INTRODUCTION 

Graphene is a new material that has been attracting a lot of attention since its experimental discovery 
[l[. This carbon aiiotrope has unique transport properties @0i hire a high electronic mobility [ij and 
thermal conductivity [5j, which are believed to be important for future applications in nano-devices 
00- However, there are certain discrepancies in the values of the electronic mobilities depending on 
whether the samples are suspended or in a substrate [HQ. At low temperatures, impurity scattering can 
be responsible for this effect, which eventually leads to a metal-insulator transition since the mobility 
edge appears near the Fermi energy [Icj . as has been confirmed in graphene doped with H (Tlj . However, 
above T > 10° K such discrepancies are believed to be a consequence of the crucial role of flexural phonon 
modes in the electron scattering, as has been shown very recently by applying tension to graphene sheets 
Q. In fact, flexural modes are collective atomic displacements which are perpendicular to the graphene's 
plane [12], making it behave like a corrugated surface. From a microscopic point of view, the scattering 
results from changes in the distances between atoms, leading to fluctuations of 7r— orbitals electron 
wavefuctions overlaps [l3| . 

However, since charge carriers in flat graphene are described by massless Dirac fermions [H, fl3 ] , it 
is natural to ask if it is possible to modify the Dirac equation taking into account the flexural mode 
interaction. Two paths can be followed in order to answer this question. One is to start from the 
usual tig ht-binding approach and use Taylor expansion of the overlap integral on the displacement 
field |l2|[l3j [l5| . Here we present an alternative point of view, in which the interaction is included by 
making the observation that a graphene membrane can be considered as a curved space. The effective 
equation must be covariant due to simple and general physical arguments. The desired equation can 
be considered as akin to the Dirac equation in curved space-time. In fact, a similar approach has been 
taken in references [l6|[l7|(l8j to study smooth ripples in graphene, using a metric that was considered 
asymptotically flat and with cylindrical symmetry. Our approach does not contain such requirements 
since the metric is general and thus can include a spatially non-decaying phonon field. This approach 
leads to interesting predictions, like the possibility of building gates for graphene electronic devices by 
superposing sinusoidal strain fields in different directions. Also, since the planar vibrational modes can 
be described by vector potential [l2| [l]| [l5[ , our approach allows to include all phonon branches and the 
electromagnetic potential in a single equation. Note that in fact, the problem of the two dimensional 
Dirac equation including a vector potential has been solved recently [19|(20(. As we shall show in the 
conclusion, the present approach has certain advantages over the Taylor expansion of the tight-binding 
parameters. 

To finish this introduction, let us briefly sketch the ideas behind the present approach. For graphene 
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at low temperatures, its surface can be considered as flat. The corresponding unperturbed Hamiltonian 
operator used to describe charges evolving on the graphene sheet can be written as an effective Dirac 
equation (j.[T^|. 

H = -thv F [ 1 x w x + 1 yw y ], (i) 

where Vf is the Fermi velocity, with 

t=-.= (?;)•' =-•=(?»)■ v 

Note that at this stage there is no difference between covariant and contravariant indices, because we 
have 

7V + 7 J 7 l = <*i<rj + o-jOi = 2g i:j 1, = x, y), (3) 

with gij = Sij , so that obviously the contravariant metric raising the indices is also g tJ = S lJ . But this is 
no more true when the underlying two-dimensional space is not flat, but corrugated, with a non trivially 
deformed metric, as the case is for graphene at T >> 10°-ftT. This is why a more careful treatment of 
this two-dimensional version of the Dirac equation should be considered. Luckily enough, the problem 
of covariant formulation of Dirac's equation in a curved space has been quite deeply investigated since 
a long time [U [22| , so that we can follow the same steps in this particular case: introduce the non- 
Euclidean two-dimensional metric on the graphene sheet, then adapt the Clifford algebra and find the 
Christoffel connection, and finally assemble all these in the covariant version of Dirac's equation. 

The layout of the present article is the following: Section II describes the induced metric due to a 
corrugated surface, and in Section III we obtain the resulting covariant Dirac equation. Section IV deals 
with the application of the generalized Dirac equation to the description of electrons evolving on the 
corrugated graphene sheet, and finally, conclusions are given in Section V. 



II. THE INDUCED METRIC 

At any finite temperature, phonons produce a displacement field (it). We start in this section by 
finding the resulting metrics and its Christoffel connection on corrugated graphene. Let the surface 
materialized by the graphene sheet be described by, 

z = f(x,y,t), (4) 

where x,y are the coordinates in the graphene sheet, and z is the out of plane displacement. Notice 
that here we only consider the flexural phonon branch, since the planar phonons can be described by 
inclusion in a vectorial potential 12], whose solution is basically known [f3|[2Cl|. 
As now the differential dz becomes a linear combination of dx and dy, we have 

dz = ^j-dx + y~dy, (5) 

and the induced metric on the sheet is given by the following formula: 

ds 2 = gij dx % dx\ = x, y, z). (6) 

where is the metric tensor. Or, in a more explicit form, 

ds 2 = dx 2 + dy 2 + dz 2 = dx 2 + dy 2 + (^J-dx + ^-dyf . (7) 

ox ay 
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and it is obviously a metric in a 2-dimensional (but curved) space parametrized two space variables (x, y) 
After opening the last expression we get the explicit form of the induced metric, which is 



~ | v dx l dx dy \ / o \ 

dx dy 




This can be written as 

9jk = 9jk + h jk , with j, k = x,y. (9) 

where gty = diag(l, 1) is the flat metric in 2-dimensional space, and h%j is the perturbation (supposed 
small as compared with 1) provoked by the corrugation of the sheet. The matrix ([8|) is symmetric and 
real, therefore it can be diagonalized by an appropriate linear transformation. 

Also the inverse (contravariant) metric can be easily found. The determinant of the matrix corre- 
sponding to the covariant metric tensor © is easily found to be 

det [ sit i ) =( - 1 + ^ / ) 2 + ^y f ^ ^ 

\ dx dy 1 + ^ dy > ) 

and the inverse matrix, corresponding to the contravariant metric g^ k is 

Jk — I ^dy> Ox dy \ f-,-i\ 

where we used the abbreviate notation for the determinant, Q = 1 + (d x f) 2 + (d y f) 2 . 

The deformation hjk is entirely composed of quadratic terms containing products of spatial partial 
derivatives of the deformation function f(x,y,t), 

h Jk ~d 3 fd k f, (12) 

and it is easy to check that the corresponding Christoffcl symbols reduce to 

f ) k = \ t m (dfgmk + d k ~g jm ~ d m g jk ^ = g m d m f d%f. (13) 

Both quantities disappear when / — 0, and the metric becomes flat again. 

The Dirac equation in two dimensions should be now generalized in order to incorporate the fact that 
the metric on the surface of constraint is no more flat, but curved. In the following section, we obtain 
the covariant generalization of the Dirac equation. 



III. THE COVARIANT DIRAC EQUATION 

The equation we want to produce now can be written as the new deformed quantum-mechanical 
Hamiltonian acting on a two-component spinor ^> as follows: 



H * 



fV 3 



+ f Vj 



(14) 



Here not only the contravariant metric is deformed, but also the 7-matrices should be modified in order 
to satisfy new anti-commutation relations with the induced metric instead of the fiat one as before; 
finally, Vj contains not only the electromagnetic and on plane phonon interaction visualized by the 
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vector potential included in the usual gauge-invariant way, but also the Christoffel symbols of the metric 

9if 

V, * - (dj - e A,-)* + f% g ki S mfc (15) 
where the Christoffel symbols T^, are defined as usual, by means of the modified metric: 

Fj-k = \ 9 im [djg m k + d k g jm - d m g jk ], (16) 
and £ mJ fc is the matrix-valued anti-symmetric tensor defined by means of the modified gamma-matrices: 

S m fe = ^[7m7t - 7fe7m]- (17) 
o 

This term is often called "spinorial connection" . [2l], [22| ■ 

We are looking for two "deformed" generators of the Clifford algebra % and 7 y that would satisfy 

fdf\\ _ dfdf . . fdf\ 2 

In order to do this, we must introduce the third Pauli matrix, because the deformation of the sheet 
pushes it out of the strict two-dimensional plane (x,y). 

The undeformed Clifford algebra satisfying anti-commutation relations in a three-dimensional flat 
space is defined as follows: 

Ijlk + lkl = 1g k 1 with g ok = diag(l, 1, 1), j, k = 1,2, 3, (19) 

is easily found to be generated by three Pauli matrices as follows: 

7i=<r x , l2 = cr y , 73 = cr z , (20) 

with 



(21) 



on _ (0 -A _ (l o 

1 oj ' ay - [i J ■' (Tz ~ \Q -1 
we have indeed 

( 7 ,) 2 = 1, (7«) 2 = 1, (7,) 2 = 1 

and the three matrices anticommuting with each other. 

The ansatz for the two deformed space-like 7-matrices is simple: if we set 

% = <?x + acr z , % = <7 y + ba z , (22) 

then the coefficients a and b should be, as it was easy to check, 

fl= fe' b= ^ (23) 

so that 

1x = cj x + — o z , 7y = °y+Q- (24) 

Now we have to produce their contravariant counterparts that appear in the Dirac equation 22]. We 
have: 

r=9 xx % + ~9 XV %, f = 9 yx % + 9 yv %, (25) 
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which gives explicitly 



T = 



Q 



1 + (dyf) 2 a x + (8 x f) a z - (d x f)(d v f) a y 



1 Q 



1 + (d x ff )a y + (d v f) a z - (d x f)(d y .f) a x 



(26) 



Recalling that Q = 1 + (d x f) 2 + (d y f) 2 , we can add and substract in the numerators of the above formula 
respectively the following terms: (d x f) 2 a x in the first one, and (d y f) 2 a y in the second one; this will 
enable us to separate the undeformed matrices a x and a y and the genuine deformation terms containing 
spatial derivatives of /. This gives the following result: 

k ~ a x V x + ay V y + ^ (grad/ ■ V) - ^ (grad/ • <?) (grad/ • V) (27) 

where the vectors and their scalar products are two-dimensional, i.e. we mean: 

grad/ = [d x f, d y /], a= [a x , a y ], V = [V x , V y ], 

so that 

(grad/) • (a) = d x f a x + d y f a y , etc. 

We see that already in the numerators we have not only linear, but also quadratic terms, notwithstanding 
the presence of quadratic terms in the denominator (contained in the normalizing factor Q). If we decided 
to keep linear terms only, then the modified Hamiltonian would contain only one extra term proportional 
to the matrix <j z : 

H Un ~ Ho + SH = a x V x + a y V y + a z (grad/ • V) (28) 

Note that also the differential operator V is taken in its primary form, because the connection coefficients 
contain only quadratic expressions in derivatives of /. Observe that this equation is akin to the one 
obtained using a Taylor expansion of the overlap integral in a tight-binding approach. 

However, if we choose to keep all terms up to quadratic ones, then we must take into account also the 
Christoffcl coefficients in V. It is easy to check the following explicit form of our Christoffcl symbols; 
keeping only the second order expressions means that we can use the simplified formula in which g 1 ^ is 
replaced by g % i = <5 U . We have then: 

r II = 9 xx f, IT ' xy = ^y X = d X f d xy f, IT = 8 X f dyyf, 

^XX ~ 3 y fd xx f, = — dyfd xy f, ~ d y fd y yf, 

In covariant derivatives, these coefficients are contracted with gi k and the anti- symmetric matrices S^ m , 

ri ~km sr 
jk9 ^im, 

which, taking into account that only diagonal terms in g lk do not vanish and are equal to one, leads to 
the following result when explicited: for j = x we have 

px xx y , px yy y , yy xx y , ra yy y 
L xx9 Z-'xx "I" 1 X y9 xy t 1 xx9 ^yx + 1 X y9 ^yyi 



(i 



and for j = y we get a similar expression: 

-px xx y I pi n yy y i -py n xx y . ™ yy yi 
iyzff ^zan+lyyff ^xy + L yx 9 ^yx + L yy g L yy , 

The non- vanishing metric tensor components are equal to one, whereas the E-matrices are anti-symmetric 
in their two lower indices, so what is left is only 

^xy ^xy + F xx ^yx = I 9 x f d xy f — d y f d xx f J Y> xy = — (dyfd 2 x f)Y> xy , 



1 yy ^xy ^ L yx ^yx 



d y fdij - d x fd 2 vv f\x xy = -(d x fd 2 vv m yx 



(30) 



These are the only terms remaining to be included in the covariant derivatives as follows: 

V, = (d x - eA x ) - (d y fd 2 xx f)Z xy , V v = (d y - eA y ) - (d x f d 2 yy f)Y> xy . 

In the Hamiltonian, they appear multiplied from the left by the corresponding 7-matrices, but here, 
evaluating the terms coming from the Christoffel connection coefficients, already quadratic in deformation 
/, we may keep only their undeformed version, which our case are just the two Pauli matrices a x and 
a v , so that the part of the deformed Hamiltonian keeping track of the Christoffel connection is 



a x V x + <J y V y = a x (8 X ~ eA x ) + o 1 



(dxf)(d 2 f)~(d v f) (d 2 J) 



E T 



+ a v {d y - eA y ) + a y 



(dyf)(d 2 x f)-(d x f) (d 2 f)) 



Taking into account the fact that 



(31) 



{&x) 2 = 1, (ov) 2 = 1) o- x a y = —a y a x , and a x Y> xy = - a y , 

and adding and substracting terms like (d y f)(d yy f) to the first expression and [d x f)(d xx f) to the second, 
we get the following invariant form of the extra terms induced by the Christoffel connection: 

(dxf)(d 2 X yf) + (d y f)(d 2 yyf) ~ (^/)(^,/) ~ {dyf) ^f) 



— a 
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(dyf) (d 2 f) + (d x fMj) ~ (dxf)(d 2 J) - (d x f)(d 2 f) 



This in turn can be written in a more compact (and elegant !) way as follows: 

~ a ■ gmd [(d x f) 2 + (dyf) 2 ] --[<?• grad/] Af 



with 



Af = d 2 xx f + d 2 yy f; 



(32) 



(33) 



It is worthwhile to note that the expression Q33p vanishes when / is a pure monochromatic wave, 
/ = Acos(ujt — K x x — K y y), but is different from zero as soon as there is a superposition of such 
expressions, e.g. for a standing wave. It is wortwhile to mention that this simple fact provides a nice 
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prediction of the present approach, since it means that one can control the electronic conductivity 
by mixing sinusoidal strain fields. For example, if two sinusoidal strains are applied in perpendicular 
directions using piezoelectric nano-devices, this will result in a very different electronic conductivity 
when compared with the case of only one sinusoidal deformation, due to the presence of the expression 
(155)) . The reason of such effect can be understood by observing that a pure sinusoidal deformation does 
not change distances inside the sheet, as can be easily verified by applying strain to a sheet of paper. If 
the strain is applied in only one direction, the sheet can be gently bended, while a second strain in a 
perpendicular direction produces real corrugation. 

Now we are able to write down the full Hamiltonian for an electron on a sheet, taking into account 
that sheet's proper motions described by the deformation from horizontal plane given by z = f(x,y,t), 
up to the second order (quadratic terms in derivatives of /): 

H~H + ^<j z (grad/) ■ (V) - i (grid/ • 8)(gvldf ■ V) 



l -a-grld[{d x ff + {d y ff}-\ 



8 ■ grad/ 



A/. 



(34) 



The normalizing factor 1/Q in front of two first contributions can be set to 1, because it contains the 
squares of derivatives of /, and if developed, will create terms of order 3 and 4 when multiplied by the 
terms behind. 

The quantum-mechanical Hamiltonian is obtained from this expression by multiplying it by —ih and 
vf- Let us define the operator of generalized momentum, 

7? = p — eA = —ih grad — eA 



where A is a vector potential that describes an electromagnetic field [19| [2(| or in plane longitudinal and 
transversal phonons 15[ . Then we can write: 



,0 7T X — ITTy 

ir x + ITTy 



v F 



(d x f)n x + (dyf)TTy 






■(d X f) TT X ~ {dyf)%y 



ihvf 



(grad/ • 8) (grad/ ■ V) + ± 8- grad [(d x f) 2 + [d y ff] - \ 



8 ■ grad/ 



A/ 



(35) 



This last equation is the generalized Dirac equation. For graphene, it is interesting to observe that in 

the Hamiltonian derived here, the terms a ■ grad [(d x f) 2 + (d y f) 2 ] /8— 8 ■ grad/ A//4 are akin to 

the ones obtained from a vectorial potential approach [l^ . However, Eq. (|35|) contains more terms, 
including a dominant linear correction in /. 



IV. THE GENERALIZED DIRAC EQUATION ON CORRUGATED GRAPHENE 

To show how Eq. (|35|) is used in the graphene case, we consider that the displacement field can be 
written using a simple set of basis functions provided by standing waves [23l | , 



/J,q>0 



^ = -^E ^(t) [Qi%cos (l?-R-^(l?)t) + sin (f-l-^it) 



(36) 



where e^(q) is the polarization vector for a wave- vector /i is the phononic branch, w M ("^) the 
dispersion relationship and = (x, y) is the position. Notice that the notation q > indicates that 
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q x > and q y > 0. Q and Q are operators given in terms of the creation (annihilation) phonon 
operators [23] ( a l^), 



*) = 1 ( (a) 



where a runs over c and s and M is the carbon atom's mass. 
Thus, for flexural phonons the function f(x, y, t) is given by, 

2 



q>0 



(c) 



COS ( 



Q^-j sin < 



(37) 



(38) 



where fi — F means that we are dealing with flexural phonons, and the phase <pq is defined as q = 
'^■K—Lj fi (~q^)t. For the flexural branch, ojF(~cf) — otF ||~^|| 2 where ctF ~ 4.6 x 10~ 7 m 2 /s. Notice that 
Qp-fi and Qp-j have length units, so f(x, y, t) has the same units, while d x f and d y f are dimensionless. 
The generalized Dirac equation only needs as input the partial derivatives of f(x,y,t), 

dxf = -j=52<lx \Q { f.i} cos ^q + Q 



q>0 



(39) 
(40) 



q>0 



For example, if we keep only the linear correction in Eq. Q35p. the Hamiltonian reads as follows, 



H = v F 



7T X — VKy 





2vf_ r ( C 



Q^sin, 



q>0 



V -Qa, TT X - q y TT y 

Using the fact that the Fermi velocity is much higher than the flexural modes' speeds, we can drop 
the time dependence and apply perturbation theory to solve the equation, as will be shown in the 
forthcoming article. 

The correction terms that results from d x f and d y f are of the order z max /A, where z max /X is the max- 
imal vertical displacement of flexural mode, and A is the wavelength of a typical ripple, usually around 
several atomic spacings ao. Thus, the energy of flexural modes correction behaves as ~ VFPF(zmax/a>o) 
in the linear approximation, the subsequent non-linear terms being proportional to the corresponding 
powers of the adimensional parameter z max /ao, as obtained in other approaches 




V. CONCLUSIONS 



We have deduced a generalized Dirac equation that describes charge carriers moving in corrugated 
graphene. To produce such equation, an appropriate metric was found and the principle of covariance 
has been applied. The resulting equation contains a linear correction plus several non-linear terms. 
Some of these non-linear terms correspond to the Taylor expansion of the overlap integral approach 12], 
while others, including the linear correction, are new. Such terms could lead to interesting effects, like 
resonances between various phonon modes, or between flexural phonons and the oscillating external 
electromagnetic field. We expect these non-linear terms to become important at higher temperatures, 
since the dependence of the amplitude of f(x, y, t) at higher temperatures contains the square root of T, 
as can be easily shown by using the equipartition of energy principle. Also, we proposed construction of 
simple electronic gate by observing that a pure sinusoidal strain field has a very different effect on the 
electron dynamics than a superposition of sinusoidals propagating in different directions. 
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